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ABSTRACT 

This  paper  reports  real  time  explosive  gas  sensing  (DNT)  in 
atmospheric  pressure  utilizing  the  noise  squeezing  effect  that 
occurs  prior  to  a  bifurcation  event.  A  noise-squeezing  controller 
based  on  the  statistics  of  phase  noise  is  implemented  using  high 
speed  Lab  VIEW™  field  programmable  gated  array  (FPGA).  A 
high  frequency  TNT-molecularly  imprinted  fixed-fixed  microbeam 
sensor  utilizes  this  nontraditional  sensing  strategy  and  performs 
DNT  sensing  at  various  concentrations.  Experiments  are 
conducted  using  both  noise-based  and  sweep-based  bifurcation 
tracking  for  a  direct  comparison.  Results  demonstrate  noise-based 
bifurcation  tracking  is  not  only  capable  of  performing  reliable 
frequency  tracking,  but  also  shows  the  method  is  superior  to  the 
bifurcation  sweep-based  tracking.  Over  three  orders  of  magnitude 
improvement  in  acquisition  rate  is  achieved,  and,  as  a  result, 
confidence  and  precision  on  bifurcation  frequency  estimation  is 
significantly  improved  over  the  bifurcation  sweep  tracking  method, 
enabling  DNT  sensing  at  concentrations  much  below  sub-ppb 
(parts-per  billion)  level. 

INTRODUCTION 

High  sensitivity  and  selectivity  and  low  cost  are  key 
considerations  for  sensors  designed  for  trace  explosive  detection 
[1].  Low  mass,  high  frequency  and  low  cost  micro/nano  sensors 
utilizing  mass  loading  of  microcantilevers  have  drawn  increasing 
attention  in  the  area  of  mass  sensing  in  recent  years  [2,3]. 
Molecular  imprinted  polymers  (MIPs)  have  become  an  attractive 
thin-film  coating  for  many  MEMS  sensors  and  a  deeper 
understanding  of  binding  sites  in  MIPs  has  been  achieved  [4,  5]. 
In  this  work,  15  nm  films  of  sol-gel-derived  xerogels  molecularly 
imprinted  for  tri-nitro  toluene  (TNT)  have  demonstrated  selectivity 
and  stability  in  combination  with  a  fixed-fixed  beam  MEMS 
sensor  [6-7].  The  sensor  was  characterized  by  parametric 
bifurcation  sweep-based  tracking  [6].  Traditionally,  mass  sensing 
using  MEMS  has  been  achieved  based  on  the  natural  frequency 
shift  due  to  an  increase  of  resonator  mass.  However,  noise  has  set 
the  limit  of  detection  for  linear  sensing  [8].  The  ability  to  track  the 
minimum  shift  of  the  natural  frequency  is  determined  by  both  the 
intrinsic  and  extrinsic  noise  of  the  system.  Dynamics  of 
parametrically  excited  oscillators  and  their  applications  have  been 
studied  extensively  [9-11]  as  well  as  successful  attempts  to 
improve  the  effective  quality  factor  of  microcantilever  arrays 
operating  in  the  linear  regime  utilizing  parametric  amplification 
have  been  made  [12,13,15].  However,  bifurcation  mass  sensing 
has  demonstrated  superior  sensitivity  in  the  presence  of 
measurement  noise  when  compared  to  linear  sensing  in  air  [14]. 
As  shown  in  figure  1,  the  critical  location  of  the  bifurcation  mass 
tracking  is  recorded  by  repeatedly  sweeping  the  frequency  towards 
the  critical  point  until  large  amplitude  occurs.  A  system  reset  is 
required  after  each  sweep  cycle  due  to  hysteresis.  This  sensing 
method  offers  higher  sensitivity  and  noise  resistivity  over  the 
harmonic  resonant  tracking  [14],  however,  the  long  settling  time 
and  its  high  dependency  on  sweep  rate  to  noise  ratio  [16-17]  make 
it  less  desirable.  In  this  work,  we  investigate  a  different  sensing 


approach  based  on  noise  squeezing  effect  in  parametric  systems 
[15,  18-20]  that  occurs  prior  to  the  critical  point.  A  Lab  VIEW™ 
FPGA  controller  is  implemented  to  keep  the  device  operating  close 
to  the  edge  of  instability  and  inhibiting  large  amplitude  grows  (Fig. 
1).  Asa  result,  close  to  three  order  of  magnitude  of  improvement 
in  acquisition  rate  is  achieved. 


Figure  1:  Schematic  of  both  noise  squeezing-based  tracking  (left) 
and  bifurcation  sweep-based  tracking  methods  (right).  The 
bifurcation  sweep  based  tracking  method  tracks  the  bifurcation 
locations  by  repeatedly  performing  frequency  sweeps  towards 
critical  point  until  large  amplitude  results.  Then  the  device  is 
relaxed  to  a  zero  stable  state  before  the  next  sweep  starts. 
However,  in  the  noise  squeezing  control  tracking  method,  the 
device  approaches  the  critical  point  until  the  noise  squeezes  below 
some  threshold,  then  feedback  control  keeps  the  device  close  to  the 
edge  of  instability  while  maintaining  small  a  response  amplitude. 

THEORY 

In  this  work,  parametric  resonance  is  achieved  by  applying 
periodic  axial  forcing  to  a  fixed-fixed  microbeam  through  an 
external  shear  piezo,  and  the  resulting  dynamic  is  governed  by  the 
Matehieu-Hill  equation  with  noise  term,  and  the  analysis  can  be 
found  in  [21].  As  shown  in  figure  2,  far  away  from  the  critical 
point  where  zero  solution  is  stable,  the  phase  appears  to  be  random 
ranging  from  -n  to  n  due  to  thermomechanical  noise.  As  it  gets 
closer  to  the  critical  point,  phase  noise  squeezes  onto  the  slow 
manifold  that  corresponds  to  the  eigenvector  direction  associated 
with  eigenvalue  with  the  smallest  magnitude.  If  it  keeps  going,  the 
phase  noise  squeezes  more  and  zero  solution  becomes  unstable, 
large  amplitude  results.  The  large  phase  variance  where  phase 
starts  to  squeeze  makes  it  a  perfect  location  for  a  feedback 
controller  operation  to  stay  close  to  the  edge  of  bifurcation  (Figure 
i). 
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Figure  2:  Transient  Response  of  a  parametrically  excited  gyro 
before  and  after  the  critical  bifurcation  point.  A  coordinate 
transformation  of  (qh  q2)  into  a  frame  rotating  at  half  the  drive 
frequency  is  used.  The  system  undergoes  a  super critical/sub  critical 
pitchfork  bifurcation.  Far  away  from  the  critical  point,  the  phase 
(in  red)  alternates  randomly  with  a  large  variance.  As  the  device 
approaches  the  critical  point,  indicated  by  the  transition  from  red 
to  green,  the  dynamics  collapse  onto  a  one  dimensional  slow 
manifold,  where  the  phase  variance  drops  dramatically,  making  it 
an  ideal  location  for  control.  As  a  consequence,  the  phase  noise 
correlates  and  the  amplitude  starts  to  grow.  Escape  to  a  large 
amplitude  (in  purple)  occurs  as  the  parameters  cross  the  instability 
boundary. 

DEVICE  AND  EXPERIMENTAL  SETUP 

The  device  under  testing  is  comprised  of  a  fixed-fixed 
microbeam  with  a  natural  frequency  of  49.33  kHz,  as  shown  figure 
3.  The  sensor  was  spin-coated  with  xerogel-based  molecularly 
imprinted  polymers  (MIPs)  that  is  highly  selective  to  TNT. 
However,  2-4  di-nitro  toluene  (DNT)  is  used  for  mass  sensing 
experiment,  the  results  can  be  inferred  to  the  sensing  of  TNT  as 
they  are  structurally  similar. 
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Figure  4:  Experimental  setup  and  FPGA  controller  design 
schematics.  The  device  is  mounted  onto  a  shear  piezo  and  is  driven 
at  ~  twice  its  natural  frequency.  The  sensor  velocity  response 
coupled  by  a  laser  vibrometer  through  an  optical  microscope.  The 
vibrometer  signal  and  the  reference  square  wave  at  half  the  drive 
frequency  are  fed  into  a  phase  lock  amplifier  (PLA).  The  FPGA 
samples  the  outputs  of  the  PLA  and  provides  feedback  control  to 
hold  response  close  to  the  edge  of  bifurcation.  The  2-4  DNT/N2  gas 
test  experiment  is  conducted  in  a  closed  chamber  at  atmospheric 
pressure  and  room  temperature.  DNT/N2  is  generated  by  heating 
the  2-4  DNT  permeation  tube  in  a  permeation  oven  at  constant 
temperature.  DNT  vapor  is  carried  by  N2  and  the  mixture  is  fed 
into  the  test  chamber  at  constant  flow  rate  controlled  by  mass  flow 
controllers.  When  the  device  is  not  under  active  DNT  test,  only  N2 
is  used. 


Figure  3:  SEM  image  of  a  450  pm  x  20  pm  x  2  pm  sensor  coated 
with  MIPs.  The  device  used  in  this  paper  is  600  pm  x  20  pm  x  2  pm. 

Shown  in  Figure  4,  experiment  is  conducted  at  atmospheric 
pressure  and  room  temperature.  The  sensor  was  mounted  on  a  shear 
piezo,  driven  by  the  function  generator  at  nearly  twice  the  resonant 
frequency  (to  drive  parametric  resonance)  at  a  fixed  voltage  of  32 
Volts.  DNT  vapor  was  generated  in  a  permeation  oven  by  heating 
2-4  DNT  permeation  tube  at  constant  temperature  mixed  with  N2. 
The  flow  rates  were  controlled  by  mass  flow  controllers  (MFC). 
The  concentration  of  the  DNT/N2  mixture  in  parts  per  billion  (ppb) 
can  be  calculated  based  on  the  flow  rate  and  oven  temperature  for  a 
specific  permeation  tube.  The  sensor  response  is  detected  by  an 
optical  laser  vibrometer  and  the  signal  is  analyzed  in  the  phase  lock 
amplifier,  which  outputs  the  amplitude  and  phase  in  reference  to  the 
reference  signal  at  half  the  drive  frequency.  The  FPGA  controller 
samples  the  phase  data  and  is  dedicated  to  the  calculation  of  phase 
variance.  Based  on  the  statistics  of  the  phase  variance,  a  feedback 
control  is  implemented  to  adjust  the  corresponding  frequency 
change  to  keep  the  device  close  to  the  edge  of  instability.  Any 
change  of  frequency  due  to  mass  loading  is  followed  by  the  noise 
squeezing  controller  immediately. 


RESULTS 

DNT  gas  sensing  experiment  was  performed  with  the  noise¬ 
squeezing  controller  at  a  low  concentration  of  0.93  ppb  DNT/N2 
mixture.  In  this  experiment,  nine  repeated  tests  of  15  minutes  of 
pure  N2  followed  by  10  minutes  of  DNT/N2  were  conducted,  as 
shown  in  Figure  5. 


Figure  5:  Figure  shows  a  constant  concentration  of  0.93  ppb 
DNT/N2  gas  experiment.  Pure  N2  was  first  introduced  for  15 
minutes  and  was  followed  10  minutes  of  DNT/N2.  The  parametric 
drive  frequency  was  lowered  by  67  ±  3  Hz  due  to  the  polymer 
absorption  of  DNT.  The  same  process  was  repeated  9  times.  The 
linear  frequency  drift  of  0.05  Hz/min  was  accounted  for  in  the  post 
processing. 

Higher  concentration  DNT  gas  experiments  were  also 
conducted.  These  experiments  were  carried  out  by  the  noise 
squeezing  bifurcation  sensing  method  and  the  bifurcation  sweep 
tracking  method  described  in  [7]  for  comparison.  Reversibility 
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was  not  able  to  obtained  due  to  the  stiction  of  DNT  at  high 
concentrations,  hence,  instead  of  introducing  N2  purging  cycle 
after  each  DNT  testing,  higher  concentration  DNT  experiment 
were  conducted  after  the  previous  concentration  reached  its 
saturation  level.  Calibration  curve  of  frequency  shift  as  a  function 
of  the  concentration  in  parts  per  billion  (ppb)  shown  in  Figure  8  is 
used  to  find  the  sensitivity  of  the  MIPs.  Frequency  stability  is 
quantified  by  Allan  Variance  [22]  of  data  highlighted  in  red  in 
Figure  6  and  Figure  7.  The  sensitivity  of  frequency  shift  to 
concentration  is  characterized  by  the  slope  of  the  calibration  curve 
in  Figure  8.  A  summary  of  the  comparison  of  the  two  tracking 
method  is  shown  in  table  1 . 


Figure  6:  DNT  gas  sensing  using  a  noise  squeezing  controller,  (a) 
shows  the  gas  experiment  with  lowest  concentration  of  1.38  ppb 
and  highest  concentration  of  18.13  ppb.  (b)  is  a  zoom  in  figure  of 
the  portion  inside  the  square  window  in  (a).  The  highlighted 
portion  corresponds  to  pure  N2  purging  after  the  1.38  ppb  DNT/N 2 
gas  test.  However,  absorbed  DNT  did  not  completely  come  off 
and  the  frequency  did  not  fully  recover  initial  starting  frequency. 
Hence,  this  explains  why  the  highlighted  data  in  red  does  not  start 
at  zero,  (b)  is  used  to  make  direct  comparison  with  the  data 
collected  from  the  bifurcation  sweep  method  (Figure  7)  at  the  same 
concentration. 


conducted  after  the  steady  state  of  lower  concentration  was 
reached. 
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Figure  8:  Calibration  curve  for  both  the  noise  squeezing  tracking 
and  bifurcation  sweep  tracking  method.  The  slope  of  linear  curve 
fit  for  each  corresponds  to  its  sensitivity. 


Table  1:  This  table  summarizes  the  results  of  both  the  noise 
squeezing  method  and  the  bifurcation  method.  Minimum 
frequency  corresponds  to  the  Allan  variance  for  each  method  and 
sensitivity  corresponds  to  the  slopes  of  the  calibration  curve.  LOD 
is  the  ratio  between  minimum  frequency  and  the  sensitivity  of 
natural  frequency. 


Noise 

Squeezing 

Method 

Bifurcation 

Sweep 

Method 

Minimum  Frequency  (Hz) 

0.045 

5.43 

Sensitivity  (Hz/ppb) 

87.62  ±3 

94.7±13.5 

LOD  (ppb) 

0.0005 

0.06 

DISCUSSION 

Experimental  results  demonstrate  that  noise-squeezing  control 
bifurcation  tracking  is  capable  of  tracking  frequency  change  due  to 
mass  loading.  The  lowest  concentration  tested  is  0.93  ppb,  and 
good  reversibility  is  demonstrated  at  low  concentration  (Figure  5). 
The  noise  squeezing  controller  estimates  the  critical  location  at  less 
than  10  ms,  while  the  average  time  for  a  bifurcation  sweeps  take 
up  to  15-20  s.  Hence,  acquisition  rate  is  increased  by  close  to  three 
orders  of  magnitude.  The  noise  squeezing  controlled  tracking 
offers  better  frequency  stability  over  the  bifurcation  sweep  tracking 
method,  leading  to  two  orders  of  magnitude  smaller  in  lowest  order 
of  detection  (LOD).  The  DNT  response  time  is  found  to  be  around 
90  sec,  which  is  much  longer  than  the  reported  time  in  [7].  This 
can  be  attribute  to  the  degradation  of  the  polymer,  since  it  was  over 
two  years  old  when  the  experiments  were  conducted  for  this  paper. 
This  issue  is  also  reflected  in  the  sensor’s  inability  to  recover  to  its 
original  state  at  high  concentration  shown  in  Figure  6  and  Figure  7. 
However,  the  DNT  response  sensitivity  is  directly  related  to  the 
sensitivity  of  the  polymer.  Hence,  confidence  of  the  noise 
squeezing  controlled  bifurcation  tracking  ability  can  still  be 
assured. 


Figure  7:  DNT  gas  sensing  using  the  bifurcation  sweep  method. 
Pure  Ns  was  introduced  after  the  first  and  last  DNT  gas  test  (2.03 
ppb  and  10.88  ppb).  Absorption  phenomenon  persisted,  and  the 
following  higher  concentration  DNT/N 2  experiments  were 


CONCLUSION 

Experiments  successfully  demonstrate  that  the  noise  squeezing 
based  sensing  is  a  superior  sensing  strategy  over  the  sweep-based 
bifurcation  tracking  method  for  real  time  DNT  explosive  sensing  at 
atmospheric  pressure.  Three  orders  of  magnitude  improvement  in 


93 


6 


acquisition  rate  leads  to  faster,  more  confident  and  more  precise 
estimation  over  the  compared  sensing  method  [23].  Even  though 
the  sensor  response  to  the  TNT  is  slow  due  to  the  degrading  of  the 
MIPs,  it  does  not  question  the  sensitivity  of  the  noise-squeezing 
controller,  since  the  TNT/DNT  sensitivity  of  the  sensor  is  highly 
depended  on  the  sensitivity  of  the  coating.  This  issue  is  expected 
to  be  resolved  with  better  coating  material  and  is  currently 
underway. 
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